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ABSTRACT 

We present the results of the "Cosmogrid" cosmological iV-body simulation suites based on the 
concordance LCDM model. The Cosmogrid simulation was performed in a 30Mpc box with 2048'^ 
particles. The mass of each particle is 1.28 x lO^M©, which is sufhcie nt to resolve ultra - faint dwarfs. 
We found that the halo mass function shows good agreement with the iSheth fc TormenI ()1999[ ) fitting 
function down to ~ lO^M©. We have analyzed the spherically averaged density profiles of the three 
most massive halos which are of galaxy group size and contain at least 170 million particles. The 
slopes of these density profiles become shallower than —1 at the inner most radius. We also find a clear 
correlation of halo concentration with mass. The mass dependence of the concentration parameter 
cannot be expressed by a single power law, however a simple model based on the Press-Schechter theory 
gives reasonable agreement with this dependence. The spin parameter does not show a correlation 
with the halo mass. The probability distribution functions for both concentration and spin are well 
fitted by the log- normal distribution for halos with the masses larger than ~ 10® Mq. 
Subject headings: cosmology: theory — galaxies: dwarf — methods: numerical — dark matter 



1. INTRODUCTION 

According to the present standard LCDM model, the 
Universe is thought to be composed pri marily of cold 
dark matter f CDM) and dark energy ([White fc ReesI 
119781 lPeacocki[l999t ). Structure formation of the Uni- 
verse proceeds hierarchically in this model. Smaller-scale 
structures collapse first, and then merge into larger-scale 
structures. 

There is serious discrepancy between the distribution 
of subhalos in galaxy-sized halos obtained by numeri- 
cal simulations and the observed number of dwarf galax- 
ies in the Local Group (jKlvpin et al.lll999t iMoore et al.l 
[T999ah . This "missing dwarf problem" is still consid- 
ered to be one of the most serious p roblem in the CDM 
paradigm (e.g. IKroupa et al.l [20T0I) . In order to un- 
derstand the origin of this discrepancy, it is necessary 
to perform high-resolution cosmological A'^-body simula- 
tions and obtain unbiased sample of galaxy-sized halos 
with resolution high enough to obtain reliable statistics of 
subhalos since the subhalo abundance sh ows large halo- 
to-halo variations (jlshivama "eFani2009bD . 

Cosmological A^-body simulations have been widely 
used to study the nonlinear structure formation of the 
Universe and have been an important tool for a better 
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understanding of our Universe. In order to study the spa- 
tial correlation of galaxies, the first cosmological A^-body 



simulations were performed in the 1970s using approxi 
mate l y 1000 particles (e.g . [Miyoshi & Ki haral 119751: iFal 
[1971 lAaTseth et al.lll97"a lEfstathiou.l979[ r Since then 



the development of better simulation algorithms and im- 
provements in the performance of computers allow us to 
use much larger numbers of particles and have drastically 
increased the resolution of cosmological simulations. 

Today, it is not uncommon that the number of parti- 
cles exceeds 10^ in high resolution simulations. In these 
works, the size of the simulation volumes is typically 
[O(Gpc)]'^ and populations of galaxy clusters, gravita- 
tional lensi ng, and the baryon acoustic oscillation are 
studied (e.g.lEvrard et al. 2002: W ambsganss ct al . 200'l; 
iTevssier et al.ll2009HKim et al.ll2009l: iCrocce et al.ll2 010'). 
The simulation results are also used to construct mock 
halo catalogues for next generation large volume surveys. 
Others use simulations of [O(lOOMpc)]'^ volumes to study 
the internal properties of galaxy-sized dark matter ha- 
los, their formatioii, evolution, and statistical properties 



e.g.[Springel et al.ll2005l:rKlvpin et al.|[20Tol:IWhite et al.l 
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Using the results of high-resolution simulations of 
small-scale structures, we can study the fine structures of 
galactic halos, the distribution of subhalos, their struc- 
tures, and their dependence on the nature of dark mat- 
ter. This information has a strong impact on the indi- 
rect search for dark matter since gamma-ray flux by self- 
annihilation is proportional to local density if we consider 
neutralino as the candidate of dark matter. Thus, we can 
restrict the nature of dark matter using the results of 
high-resolution simulations of small-scale structures and 
indirect searches of dark matter. In addition, galaxies 
are considered to form in d ark matter halos with a mass 
large r than a critica l value (Strigari et al."2008'; 'Li et al.l 
120091: iMaccio et al.H 2009: Okamoto & Frcnk 2009). The 
structures of smallest halos which can host galaxies is 
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important for the understanding of the galaxy formation 
processes. 

The simulation of smaller-scale structures of dark mat- 
ter halos is not a trivial task since a very wide dynamic 
ranges of space, mass, and time must be covered. In 
particular, the number of time steps of such simulations 
is significantly larger than that of larger-scale simula- 
tions since the dynamical time scale is proportional to 
l.Q/^/Gp, where p is the local density. Structures of 
smaller-scales form earlier, and thus have higher densi- 
ties, therefore, simulations of smaller scales are compu- 
tationally more expensive. 

Recently, simulations with gala ctic halos of very high 
resolution have been performed 



Diemand et"aITl2o68l: 



iSpringel et all [20081: iStadel et al.l |2009D. The"se works 
used the re-simulation method, where one selects one or 
a few halos at z = from a simulation which covers a 
large volume [typically a cube of size O(lOOMpc)] with 
a relatively low-resolution. The corresponding regions of 
these halos are then identified in the initial particle dis- 
tribution, and the particles in these regions are replaced 
by a larger number of smaller particles. After this is 
done, the entire volume is simulated to z = again. 

With this re-simulation method, we can resolve 
the structu res of selected halos with extremely high 
resolution (iDiemand et al.l l2008t ISpringel et all 120081: 
IStadel et al.ll2009D . However, this method cannot be 
used for the study of halo-to-halo variations. Different 
halos are born in different environments and grow differ- 
ently. The difference in the environment and growth his- 
tory must be the cause of halo-to-halo variations. There- 
fore, in order to study variations, we need a bias-free set 
of a large number of halos. Clearly one cannot obtain 
a large number of halos with re-simulation method in 
practical time. 

In principle, one can improve the statistics by increas- 
ing the number of halos selected for re-simulations. In 
order to avoid the selection bias, we need to apply ran- 
dom, bias-free selection, and the most reliable bias-free 
selection is to select all halos, in other words, to simu- 
late the ent ire simulation b o x with uniformly high mass 
resolution. Ilshivama et al.l (|2009bl ) performed the first 
bias-free high resolution simulation of small-scale struc- 
tures. They analyzed the statistics of the subhalo abun- 
dance using the complete set of halos in the simulation 
box. The number of particles was 1600'^ in a 46.5Mpc 
cubic box and the mass of a particle was IQ^Mq. The 
subhalo a bundance showe d large halo-to-halo variations 
[see a lso Ilshivama et al.l (|2008D : iBovlan-Kolchin et al.l 
(|2010f )]. The concentration parameter and the radius at 
the moment of the maximum expansion showed fairly 
tight correlation with the subhalo abundance. Ha- 
los formed earlier have smaller number of subhalos at 
present. This correlation suggests that the difference 
in the formation history is the origi n of the v a riatio n 
of the subhalo abundance [see also iGao et all (l200l: 
Ivan den Bosch et al.l (l2005D: IZentner et al.l (I2005ll1. 

Th e Millennium-II simulation (jBovlan-Kolchin et al.l 
|2009( ) used a 137Mpc cubic box and the particle mass of 
~ 9.45 X lO^Af0. Its result is suitable for the analysis of 
the statistics of galaxy-sized dark matte r halos, because 
the num ber of halos is larger than that of llshivama et al] 
(j2009b( ). However, due to the lack of the mass resolution. 



it cannot be used to study the statistics of dwarf-galaxy- 
sized halos and the statistics of subhalos with the size 
larger than faint dwarf galaxy. 

In this paper, we describe the first result of our Cos- 
mogrid simulation. We simulated the evolution of halos 
in a 30Mpc cubic box using 2048'^ particles. The mass 
of one particle is 1.28 x IO^Mq. The resolution reaches 
down to ultra- faint dwarf-galaxy-sized halos (~ IO^Mq) 
and is more than 8 times better than t hat of our pre- 
vious simulation ("I shivama et al.l r2009b[ ). We focus on 
the halo mass function with the mass down to 10 ''Mq, 
the structures of most massive halos, and statistics of 
the internal properties of dwarf-galaxy-sized halos. We 
describe our initial conditions and numerical settings in 
Section [21 and results in Section [3l We discuss and sum- 
marize our results in Section [Jj 

2. INITIAL CONDITIONS AND NUMERICAL METHOD 

The cosmological parameters adopted are based on 
the concordance LCDM cosmological model (fio = 0.3, 
Ao = 0.7, h = 0.7, fjg = 0.8, n = 1.0). These val- 
ues a re the same as those u sed in our previous simula- 
tion (jlshivama et al.l[2009b[ ). We used a periodic cube 
of the comoving size of 30Mpc. The number of parti- 
cles for the largest run is 2048"^ which corresponds to a 
mass resolution of 1.28 x lO^M©. To generate the initial 
particle distribution s, we used the MPGRAFIC package 
(jPrunet et al.l [20081 ) . w hich is a parallelize d variation of 
the GRAFIC package (|Bertschingeill2001h . The initial 
redshift was 65. 

In order to investigate the effect of the mass and spa- 
tial resolution, we performed two simulations with lower 
resolution. We generated the initial conditions for these 
low-resolution runs (CG1024 and CG512) by replacing 
eight or 64 particles in the high-resolution initial condi- 
tion (CG2048) with a single particle eight or 64 times 
more massive. In order to study the effect of cosmo- 
logical parameters, we performed an additional simula- 
tion with cosmological parameters close to recent obser- 
vational v alues [On = 0.27, An = 0.73, h = 0.7, gs = 0. 8, 
n = 0.96, ISpergel et all (|2007| ): iKomatsu etaP (|2009f )]. 
In table [TJ we summarize parameters used in our simu- 
lations. 

We used a leapfrog integrator with shared and adap- 
tive time steps. The step size was determined as 
min(2.0-\/e/|ai|, 2.0£/|wi|) (minimum of these two values 
for all particles). The gravitational plummer softening 
length £ was 175pc at z = 0. The softening was constant 
in comoving coordinates from z = 65 (initial condition) 
to z — 10. From z = 10 to z = 0, it was constant in 
physical coordinates. This procedure is similar to that 
used in Kawai et al.! (I2004D . 

For the largest simulation, we used four supercom- 
puters. Three of them are Cray XT4 machines at the 
Center for Computational Astrophysics of National As- 
tronomical Observatory of Japan, the Edinburgh Paral- 
lel Computing Center in Edinburgh (United Kingdom) 
and IT Center for Science in Espoo (Finland). The 
fourth machine is an IBM pSeries 575 at SARA in Am- 
sterdam (the Netherlands). Part of the calculation was 
done in a "grid" computing environment, in which we 
used more than one machine simultaneously for one run 
(jPortegies Zwart et al.ll2010D . 

For the time integration we used the GreeM code 
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Table 1 

Run parameters. Here, N, L, e, and m are the total number of particles, the box length, the softening length, the mass resolution. The 
cosmological parameters are f2o = 0.3, Aq = 0.7, h = 0.7, as = 0.8, n = 1.0 for (a), and Qq = 0.27, Aq = 0.73, h = 0.7, erg = 0.8, n = 0.96 

for (b). 



Name 


N 


L(Mpc) 


£(pc) 


m(M0) 


Cosmology 


CG2048 


2048^ 


30.0 


175 


1.28 X lO'' 


a 


CG1024 


1024^ 


30.0 


350 


1.03 X 10*^ 


a 


CG512 


512^ 


30.0 


700 


8.21 X lO'^ 


a 


CG512WMAP 


512^ 


30.0 


700 


7.39 X lO'^ 


b 


IFM2009 Qshivama et al. 2009b') 


1600^ 


46.5 


700 


1.00 X lO'^ 


a 



Table 2 

Global parameters of three most massive group sized halos at z = 0. Here, M, N, i?viri flvmax, and Vmax are the mass, the number of 
particles, the virial radius in which the spherical overdensity is 101 times the critical value, the radius where the rotation velocity is 

maximum, and the maximum rotation velocity, respectively. 



Name Run j\/(lO"M0) N Rvir(kpc) Rvmax(kpc) KnaxCkm/s) 

GPl CG2048 5.24 408499843 969 200 596 

CG1024 5.19 50632942 966 186 589 

CG512 5^22 6361253 968 184 596 

GP2 CG2048 3.58 279382586 854 305 476 

CG1024 3.57 34836692 853 279 472 

CG512 3^57 4347651 852 294 475 

GP3 CG2048 2.25 175752770 731 178 434 

CG1024 2.26 22072073 732 187 431 

CG512 2.25 2746874 731 192 434 



(|Ishivania et al.l l2009al ) for single supercom puter runs 
and the SUSHI code (IGroen et al.l MTU) for multi- 
supercomputer runs. The GreeM code is a massively par- 
allel TreePM code based on th e parallel TreePM code of 
lYoshikawa fc Fukushigd (|2Q05D for large cosmological N- 
body simulations. The long ran ge forces are calculated 
by the PM method (iHocknev fc E astwood 1981), and the 
short range forces arc calculated by t he Barne s Hut Tree 
method (Barnes & Hut 198^. lYoshik awa fc Fukushigd 
(f2005i) used a 1-D slab decomposition, but in GreeM we 
use a 3-D multi-section decomposition (jMakinoll2004D to 
improve its scalability. In addition, the decomposition is 
based on CPU time measurement, so that near ideal load 
balance is archived. The SUSHI code is an extension of 
the GreeM code which can run concurrently on multiple 
superc omputers. It uses t he MP Wide communication li- 
brary (|Groen et alj I2010D to facilitate message passing 
between distributed supercomputers. We used 512'^ PM 
grid points for PM calculations, the opening angle for 
the tree method was 0.3 from initial to z = 10, and 0.5 
from z = 10 to z = 0. 

The calculation time was ~180 seconds per step with 
1024 cpu cores for the largest run on the Cray XT4 in 
Japan and ~140 seconds per step with 2048 cpu cores on 
the IBM pSeries 575 in the Netherlands. We spent about 
3.5 million CPU hours to perform all the 60283 steps in 
our simulation. 

The method of identif ying halos and s ubhalo s is the 
same as that described in llshivama et all (l2009bD. which 
is ba sed on the spherical overdensity (|Lacev fc Cold 
119941 ). The mass of the most massive halo is 5.24 x 
10^^ Mp^. It contains 4.08 x 10^ particles. Via Lactea 
I, II (|Die mand et"al l 12007 1 [200l . and Aquarius simula- 
tions (Springel et al. 200l) used ~ 10* million 5 x 10*, 
and ~ 10^ particles for the largest single halos. Table [5] 
shows the properties of the three most massive halos in 
our simulation. The virial radius of a halo is defined 



as the radius in which the spherical overdensity is A.{z) 
times the critical value. The overdensity A(z) is given 
by the analytic formula (jBrvan fc Normanlll998[) . 

A(z) = (187r2 + 82x- 39x2)/r2(z), (1) 

where x = il{z) — 1. The mass of a halo is defined as 
interior mass within the virial radius. 

Figure [T] shows the snapshots at z = 0. In Figure [H we 
also present the time evolution of the whole box and that 
of the most massive halo. The three most massive halos 
in simulations with three different resolutions are shown 
in Figure [3l The positions of subhalos agree very well 
in three simulations. Of course, there are some discrep- 
ancies near the centers of halos. In particular, whereas 
there is only one core in the center of the second massive 
halo (GP2) of CG2048, there are two cores in GP2 of 
CG1024 and CG512. 

The reason of this difference is that the formation his- 
tory of this halo is rather violent. It experienced many 
mergers near 2 = in the center of the halo and is far 
from the relaxed state. The difference of the accuracy 
of integration changed the timescale of the mergers of 
the halos with three different resolutions. At z = 0, the 
halo GP2 has just completed the merger in the CG2048 
run, whereas the same merger event is still on-going in 
CG1024 and CG512 runs. If we consider the spherically 
averaged density profile of the halo, the difference be- 
comes important (see Section [X^ . 

3. RESULTS 

3.1. Mass Function 

IPress fc Schechteil (|1974[ ) established a recipe to de- 
rive the number of dark matter halos based on the hi- 
erarchical clustering model. Since then, a number of 
analytic formulae for the mass function have been pro- 
posed. Many of them are designed to give a good 
agreement with results of high-resolution iV-body sim- 
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Whole resit>n « 




Figure 1. The projected density of dark matter at 2 = in our largest simulation (2048^ particles). Top panel shows the whole region 
with the volume of (SOMpc)'^. Bottom panels show the projected density of the two most massive group sized halos. These volumes are 
(2Mpc)3. 
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Figure 2. The evolution pictures of our largest simulation. Top six panels show the evolution of the whole region. Bottom six panels 
show the evolution of the most massive halo. 
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GPl CG2048 



GPI CGI024 



GPl CG512 




Figure 3. The projected density of dark matter at 2 = 0. Eacli row sliows one of the three most massive halos with mass decreasing 
from top to bottom. Columns show different resolution from highest (left) to lowest (right). 



ulations (^e.g. iSheth & Tormen 


1999 


Jenkins et al. 


2001 


Reed et al.' 2003: Yahari et al. 


2004 


Warren et al. 


2006 


Tinker et al..,2008t and references therein"). 



These formulae can reproduce the mass function be- 
tween 10"'^'' M0 and cluster scale very well. Here, we ex- 
amine the mass function of mass below lO^^M© down to 
IO^Mq. The mass functio n of this range has been studied 
only in high redshift (e.g. iReed et aITl2007l : iLukic et"al] 
[200% . 

Figure S] shows the halo mass functions at three dif- 
ferent rcdshifts for CG2048 ru n and the prediction o f 
Sheth & Tormen formula (ST. ISheth fc TormenI [19991 ) . 

The agreement is very good for the mass from ^ IO^Mq 
to M 1.0 X lO^M© at z = 0. The difference is less 
than 10% for M = 5.0 x 10^A/q to M = 2.0 x lO^^^fo 
at z = 0, M = 5.0 X IO^Mq to M ^ 5.0 x IO^^Mq at 
z = 3, and M = 8.0 x lO^Mm to M = 4.0 x lO^M^ at 



z = 5.4. 

Clearly, the difference becomes larger for small mass 
limit as the redshift increases. This might be caused by 
the softening length we used. The virial radius of halos 
becomes smaller as the redshift increases. In other words, 
the ratio -Rvir/e is smaller for high redshifts. Thus, the 
more redshifts increase, the more the softening may re- 
duce the number of small halos. 

Our results imply that the mass function is well repre- 
sented by the ST function down to 10 ''M©. However, our 
simulations have a slightly larger number of halos than 
the number predicted by the ST formula. Note that the 
finite volume of our simulation (the box length is SOMpc) 
might affect the mass function in some degrees. The ab- 
sence of long wavelength perturbations might increase 
the number of intermediate mass halos by about 10% 
(|Bagla fc Prasad|[2006l: |P ower fc Knebel[2006l) . 
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10^ 10' 10* 10' 10'°10"lO"lO"lO"' 10' 10' 10* lO'lO^'lO^'lO^lO^lO"* 10" 10' 10* 10' 10'° 10" 10" 10" 10" 

M [h 'M,, ] M [h 'M . ] M [h 'M,. ] 

Figure 4. The mass func tion of our largest simul ation (CG2048). The results of z = 0.0 (left), z = 3.0 (middle), and z = 5.4 (right) are 
shown. Solid curves are the lSheth fc TormenI II1999I ) function. Error bars are Poisson errors. 
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Figure 5. Spherically averaged radial density profiles of largest 
three halos at z = 0. Two of three profiles (middle and bottom) 
are vertically shifted downward by 1 and 2 dex. Vertical dashed 
lines show the softening length of three simulations. Upside short 
vertical bars indicate the reliability li mit of the most massive halo 
calc ulated using criterio n proposed bv lFukushiee fc Makinj 11200 II ) 
and I Power et at] 120031) . The red, blue, and green correspond to 
the simulation CG2048, CG1024, and CG512. 

3.2. Density Structures of Most Massive Halos 

Many groups have studied the density profile 
of dark matter halos using high-r esolution cos- 
mological A^-body simulatio ns (e.g . iNav arro ct al 
1997t iFukushige fc Makino( 119971: iMoorc ct al 

Jill, 
20011 



iGhigna et al.l 



- 200(T I.Ting fc Suto 2000 

2000t IFukushige fc Makhlol 120011: i Klvpin et al.. . _ _ . 
Taylor fc Navarrol 120011: jjhi g & Sutd 120021 iPower et al 



2003 : 
2004; 



I 



Fukushige & Makino 2003: 'Fukushige et al 
Diemand ct al. 2004; Havashi ct al. 2004, 
120041: iDiemand et all 2005; Reed et al 
'20061: iDiemand et al, ,200 



Navarro et al 
200,1 iKazantzidis et al 



Gao et al. iStadel et al. I l2009HNavarro et al.ll201(f . 

In most of recent works, the slopes of radial density 
profiles were around —1 in the inner region and around 
—3 in the outer region. The slope of density became 
shallower as the radius becomes smaller. Thus, the 
central slope is not described by any single power. 
Furthermore, the density profile was not universal. In 
other words, the slope showed a significant halo-to-halo 

scatter. 

Recent studies (jStadel et al.ll2009l: INavarro et al.ll2010l) 
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Figure 6. Slopes of radial density profiles of largest three halos at 
z = 0. Top panel shows those of the largest halo for three different 
resolution. Bottom panel shows those of largest three halos for the 
largest simulation (CG2048). 



based on high resolution simulations of galactic halos 
showed that the slopes of density became less than — 1 
at the radius 0.001 times the virial radius of the halo. 
Einasto profile showed bett er agreement than the NFW 
profile (|Navarro et al.|[r997f) which has been widely used 
for modeling dark matter halos because of its simplicity. 

Almost all recent high-resolution simulations of single 
halos used galaxy-sized halos. Therefore, little is known 
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if these finding can be applied to halos of different masses. 
Here, we present the density profiles of three most mas- 
sive halos in our simulation. These halos are galactic 
group-sized ones, with the mass of 5.24, 3.58, and 2.25 
xIO^^Mq. They contain 408, 279, and 176 million par- 
ticles. 

The top panel of Figure [5] shows the spherically aver- 
aged density profiles of these halos at z = 0. We can 
see that the results of three simulations with different 
resolution are indistinguishable for radii larger than the 
reliability limits, except for the second massive halo. We 
cal culated the reliability limits usin g criterion pr oposed 
bv iFukushige fc Makinol (|200l and lPower eFaLl ([2003. ). 
We can not ignore the effects of the local two-body relax- 
ation for radii smaller than these limits. As can be seen 
in Figure [3l the slight difference of the merging epoch of 
the central cores caused this difference. 

The slopes of density profiles become gradually shal- 
lower as the radius becomes smaller. The top panel of 
Figure [6] shows the slopes of density profiles of the most 
massive halo. As in the case of the density profile itself, 
the slopes also agree well with each other. The bottom 
panel of Figure |6] shows the slopes of the three most mas- 
sive halos in CG2048 run. 

The slope at 0.001i?vir is -0.9 1.0 for our three 

halos. This value is in excell ent agreement with the 
result of Aquarius simulation ("Springel et al.l I2008D or 
GHALO simulation (StadeLel al. 2009). Both of them 
gave the slope —1.0 for r = 0.001r2oo- this agreement 
does not mean the density profile obtained by these 
simulation and those by our simulation are identical. 
The concentration parameter, which we define here as 
Cvmax = i?vir/-Rvmax is 4.8, where i?vir and i?vmax are 
the halo virial radius and the radius of the maximum 
rotational velocity. This value is significantly smaller 
than that of Aquarius A-1 halo. Thus, the Aquarius 
halo is significantly more centrally concentrated, and yet 
the slope at r = 0.001i?vir is the same. Thus the rate 
of the shallowing of the slope is somewhat faster for the 
Aquarius halo than for our CG2048 halos. Most likely, 
this difference is due to the difference in the mass of the 
halo. 

In the Cosmogrid simulation, there are 8 halos with 
more than 100 million particles, and 54 halos with more 
than 10 million particles. We will analyze the density 
profiles of these halos to study the dependence on the 
mass and other parameters. 

3.3. Concentration Distributions 

The concentration parameter has been widely used to 
describe the internal structure of halos si nce it is tightly 
corre lated with the formation epoch (jWechsler et al.l 
|2002| ). Usually, the concentration is parameterized as- 
suming that the density profiles of halos can be fitted by 
the NFW profile (Navarro ct al. 1997), 



p{r) = 



Pa 



{r/rs)il + r/r,y 



(2) 



where po is a characteristic density and is a scale ra- 
dius. The concent ration c ^fw = R /ts is widel y used 
Ce.gjBullock et al.|[2 001b: Z hao et al.1' 2003: Macci 5 et al 
2007; 'Neto et al.' 2007; Maccio et al.l 120081: I Zhao et al 
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2009;.Mufioz-Cuartas ct al...2010i) . It is known that cnfw 



Figure 7. Residuals of concentration Cvmax = R-vir/R-vmax from 
tlie largest simulation (CG2048) to the lower resolution simulation 
(CG512). 

depends weakly on the halo mass. Halos with higher 
mass have smaller concentration, since the average den- 
sity of a halo refiects the cosmic density at its forma- 
tion time. The de pendence is weaker for higher redshift 
(jZhao et al.|[200l . 

The concentration based on the NFW profile is a ffected 
by fitting ranges and resolution (iNeto et al.ll2007t ). Fur- 
thermore, recent high resolution simulations showed that 
the den sity profile is significantly different from the NFW 
profile ()Stadel et al.l 120091: iNavarro et all 120101 also see 
Section [5^ . T hus, the use of Cnfw might cause some 
systematic bias (|Gao et al.ll2"008l: iReed et al.ll2010[) . 

We use the concentration Cvmax defined in Section [3.21 
which is a simpler quantity to measure the concentration. 
Note that i?vinax can be easily determined directly from 
spherically averaged mass distribution without the need 
of any fitting formulae. If the density profile is repre- 
sented by the NFW profile, either concentration can be 
converted to the other. 

First, we determine the minimum number of particles 
in a halo necessary to reliably determine the concentra- 
tion. Figure [7] shows the normalized difference of average 
concentration between the G2048 run and the CG512 run 
as the function of halo mass. We can see that the differ- 
ence is ^ 0.05 for the halo mass larger than 3.0 x lO^^Af©. 
For halo mass less than 3.0 x IO^^Mq, the difference is 
larger. In the CG512 run, a halo of mass 3.0 x lO^^Mg 
contains ^ 4000 particles. So we conclude that we need 
^ 4000 particles to reliably determine the concentration. 
For the CG2048 run, the rehability limit is 5.0 x lO^M©. 

Figure [8] shows the median, and first and third quan- 
tiles of the concentration as a function of the virial mass 
of the halo. We can see a clear correlation between the 
halo mass and the concentration. Apparently, the depen- 
dence is weaker for smaller mass. Therefore, the fitting 
functions with a single power (e.g. iBullock et al.ll20"01al : 
INeto et al.ll2007l : iMaccio et al.ll2007t iKlvpin et al.ll2010t ) 
cannot be used for halos of the size of dwarf galaxies. 

Theoretically, the concentration of a halo refiects 
the cosmic density at the formation time of the halo 
(|Bullock et al.l l2001aD . The concentrations of halos 
formed earlier are higher than that of halos formed later. 
However, the dependence should be weak for small halos 
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Figure 8. Concentration plotted against tlie lialo virial mass M 
at 2 = 0. Circles show the median value on each bin. Whiskers 
are the first and third quantiles. The number of halos on each bin 
are shown below circles. Thick solid line shows the result from an 
analytical model. 

since the dependence of the formation epoch to the halo 
mass is small for small (smaller than IO^AIq) halos. The 
slope of the power spectrum of initial density fluctuations 
approaches to —3 for small mass limit. 

In Figure [HI we also plot an analytical prediction of 
the ma ss-concentration relat ion, obtained by the method 
used in lNavarro et al.l (|1997D assuming that all halos have 
the NFW profile. The formation redshift Zf of halos with 
the mass M is defined as the epoch at which progenitors 
with the mass larger than fM first contained the half of 
the mass M. It is estimated by usin g the Press Schechter 



py • 

formalism (e.g. lLacev fc Colelll993| ) 



erfc • 



1 

2' 



(3) 



where Scnt{z) is the critical overdensity for the spherical 
collapse at Z{, and CTq(M) is the variance of the density 
fiuctuation at z = smoothed by a top-hat filter on a 
mass scale of M. Here, we used / = 0.01. The char- 
acteristic density po of a halo should refiect the cosmic 
density at the formation time. Thus, we assume 

Po = Pnorm (1 + ^^f)^ , (4) 

where Pnorm is chosen to fit the simulation results. The 
mass of a halo with the NFW profile is given by 

M = 4TTpor^AHi + c)-c/{l + c)]. (5) 



The mass and concentration at z 
other by 



are related to each 



4 4 
M = -7ri?3.^A(0)pcrit = ^Trrrc 



3^3 



A(0)p, 



crit ; 



(6) 



where pcrit is the critical density. From equations ([T]), 
©J (HI, (ISI), and we can analytically estimate the 
concentration of ha los with the m a ss M . 

As mentioned by iLacev &: Cold (|1993D . the estimated 
formation epoch obtained using equation ^ is not nec- 
essarily correct. Nevertheless, as seen in Figure |51 the 
analytical prediction based on equation ([3]) shows a very 
good agreement with the result from CG2048 run for 
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Figure 9. The probability distribution functions of the concen- 
tration at 2 = 0. These panels show the results of different mass 
ranges. Dashed curves are the best fits of the log-normal distribu- 
tion. 

halos with mass smaller than 10 ^^M0. For halos with 
the mass larger than IO^^Mq, the difference between 
CG2048 results and analytical ones are relatively large. 
However, these halos are rare objects in CG2048 run, 
and the fact might affect the results in some degrees. 
We can conclude that the shallowing slope of the mass- 
concentration relation naturally emerges from the nature 
of the power spectrum of initial density fluctuations. 

The slope is slightly shallower than that of cnfw 
for larger halos. For the case of Cnfw, the slope is 
aro und —0.10 for re l axed halos and —0 .11 for all ha- 
los ()Neto et all 120071: iMaccio et aI1l2007[ ). On the other 
hand, for the CG2048 simulation, the slope is around 
-0.07 for halos with the mass IO^^Mq, and -0.06 for 
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Figure 10. Residuals of spin from tfie largest simulation 
(CG2048) to lower resolution simulation (CG512). 
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Figure 11. Spin parameter A plotted against the halo virial mass 
M at z = 0. Circles show the median value on each bin. Whiskers 
are the first and third quantiles. The number of halos on each bin 
are shown below circles. 

halos with the mass lO^M©. Note that one overesti- 
mates the central density of halos if one estimates the 
concentration of dwarf-sized halos by extrapolating the 
mass-concentration relation of galaxy or cluster-sized ha- 
los. 

Figure [5] shows the probability distribution functions 
of the concentration parameter at z = in two different 
mass ranges. Both shapes are well fitted by the log- 
normal distributions, 



1 



2tt<t 



exp 



log^ (Cymax/Cp) 

2a2 



(7) 



We find logco = 1.050, cr = 0.124 for halos with the 
mass of 5.0 x IG^Mq <M< IO^Mq, logco = 1.022, a = 
0.128 for halos with the mass of IO^Mq <M< 10^°Mq, 
and log Co = 0.965, a = 0.125 for halos with the mass of 
10i°A/(^ < M < 10"Mp 



/0. 



3.4. 



Spin Distributions 

The dimensionless spin parameter is a good parameter 
to quantify the rotation of a halo. One often uses the 



spin parameter defined in iBullock et al.l (|2001aD , 



^/2MVR' 



(8) 



where M, R, V, and J is the virial mass of the halo, 
radius, rotational velocity at R, and total angular mo- 
mentum inside R. 

The distribution, the dependence on the halo 
mass, and the evolution ha ve b een studi ed by 
a number of works (e.g. [Bul lock ct ahl I2001al: 



Bailin & Steinmetz T005; Bett et al. 2007 : iMaccio et alj 

2007; Kncbc & Power 

Antqnuccio-Delogu et al.l 
20101: IWang et al.l I2010D 



Macc io et al.l 120081: 



2008; 

120101: iMunoz-Cuartas et al.l 
The spin of galaxy sized 
halos are well studied by using the results of sufficient 
resolution simulations. However, we do not understand 
those of dwarf galaxy sized halos. The spin distribution 
of those halos at only high redshifts ar e studied by the 
result of high resolution simulation (|Knebe &: Poweii 
[200l) . Here, we extend the spin distributions at z = 
to dwarf galaxy sized halos (down to IO^Mq) in the 
same way as the concentration. 

First, we determine the minimum number of particles 
in a halo necessary to reliably determine the spin as done 
for the concentration. Figure [TU] shows the normalized 
difference of average spin between the CG2048 run and 
the CG512 run as a function of halo mass. We can see 
that the difference is ^ 0.05 for halo mass larger than 
8.0 X IO^Mq. For halo mass less than 8.0 x IO^Mq, the 
difference is large. In the CG512 run, a halo of mass 
8.0 X 1O^M0 contains ^ 1000 particles. So we conclude 
that we need ~ 1000 particles to reliably determine the 
concentration. For the CG2048 run, the reliability limit 
is 1.28 X 10^ Mq. 

Figure [Tl] shows the median, and first third quantiles 
of the spin parameter as a function of the virial mass 
of the halo. Apparently, we can see the spin parameter 
is independent of the mass down to lO^M© as pointed - 
out for larger halos in previo us works ([Maccio et al.ll2007l : 
iMufioz-Cuartas et al.ll20i ^. The median value is 0.0336. 

Figure fT2l shows the probability distribution functions 
of the spin parameter at z = in two different mass 
ranges. Both distributions are almost identical. The dis- 
tributions are well fitted by the log-normal distributions. 



P(logA) 



27rCT 



exp 



log' (A/Aq) 
2a2 



(9) 



We find logAo = -1.477, = 0.308 for halos with 
the mass of 1.28 x IO^Mq < M < lO^M©, logAo = 
-1.480, cr = 0.288 for halos with the mass of lO^Afg < 
M < IO^OMq, and logAo = -1.472, cr = 0.277 for halos 
with the mass of lO^^M© < M < 10"Mq. Thus, we 
conclude that there is no mass dependence of the spin 
parameter. Otherwise, it is extremely weak. 

We can see that there are small deviations from 
the log-normal distributions at high spin regions as 
seen i n previous works for la r ger h alos ([Bett et al.l 
I2OO7I : lAntonuccio-Delogu et all mM- ^ 

we re- 

move unrelaxed halos, the de viations will disappear 
([Antonuccio-Delogu et al]|2010D . 

4. DISCUSSIONS AND SUMMARY 
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Figure 12. The probability distribution functions of the spin 
parameter at 2 = 0. These panels show the results of different 
mass ranges. Dashed curves are the best fits of the log-normal 
distribution. 

We present the first scientific results of the Cosmogrid 
simulation. Because of unprecedentedly high resolution 
and powerful statistics, the simulation is suitable to re- 
solve internal properties of halos with the mass larger 
than dwarf galaxy and subhalos whose scales are compa- 
rable to ultra-faint dwarf galaxies. 

We summarize the main results of this paper as follows: 

• The halo mass f u nction is well described by the 
iSheth fc TormenI (jl999( ) fitting function down to 
~ IO^Mq from 1.0 x IO^^Mq. The differences are 
less than 10% at z = from M = 5.0 x lO'^Af© to 
M = 2.0 X IO^^Mq. 

• We analyzed the spherically averaged density pro- 



files of the three most massive halos which contain 
more than 170 million particles. Their mass are 
5.24, 3.58, and 2.25 xIO^Mq. We confirmed that 
the slopes of density profiles of these halos become 
shallower than —1 at the inner most radius. The 
results are consistent with the recent studies based 
on high resolution simulations for galactic halos. 

• We studied internal properties of halos at z = 
with the mass more than ~ lO^M©. The concen- 
tration parameter measured by the maximum rota- 
tional velocity radius is weakly correlated with the 
halo mass. We found that the dependence of the 
concentration parameter with halo mass cannot be 
expressed by a single power law, but levels off at 
small mass. The slope of the mass-concentration 
relation is around —0.07 for halos with the mass 
lOi°Af0, and -0.06 for halos with the mass 10^ Af©. 
The shallowing slope naturally emerges from the 
nature of the power spectrum of initial density 
fiuctuations. A simple model based on the Press- 
Schechter theory gives reasonable agreement with 
the simulation result. The spin parameter does not 
show a correlation with the halo mass. The prob- 
ability distribution functions of concentration and 
spin are well fitted by the log-normal distribution 
for halos with the mass larger than ^ lO^M©. 

We have shown here a first analysis of the Cosmogrid 
data and we plan to extend our analysis in future pub- 
lications. Some of the topics that we want to address 
are: the variation of density profiles and its impact on 
the dark matter detectability, the statistics of subhalo 
abundance of the mass scale down to ultra faint dwarf 
in dwarf- and galaxy-sized halos, the assembly histories 
of halos, and the evolution of internal properties of halos 
which are presented in this paper at only z = 0. 
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